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Possibility of structured mixed phases at first order phase transitions is examined by 
CSJ . taking into account of charge screening and surface effects. Hadron-quark phase transition 



in dense neutron star interiors is considered as a concrete example. 
1. INTRODUCTION 



Recently many efforts have been made for the derivation of the equation of state (EOS) 
of the neutron star matter exhibiting possibilities of various phase transitions in neutron 
- star interiors. Some of these possible phase transitions are the first order phase transitions 
! (FOPT) essentially affecting the EOS. We will discuss this issue by considering the hadron- 
F^. ', quark (H-Q) matter phase transition [ |l]]. 

In Ref. [ Glendenning has pointed out that the usual Maxwell construction is appli- 
cable only for systems with one particle species and respectively one chemical potential, 
whereas in neutron stars there are two relevant quantities, the charge and baryon number 
chemical potentials. He has advocated to use the Gibbs conditions and demonstrated the 
appearance of a wide region of the structured mixed phase, consisting of droplets, rods or 
slabs of one phase embedded into the other phase, instead of the Maxwell construction. 
Finite size effects, the Coulomb interaction and the surface tension, in the treatment of 
the mixed phase were disregarded. His idea was further exploited by many authors. It 
was concluded that a wide region of the mixed phase is inevitable and there should be no 
region of constant pressure, opposite to that follows from the Maxwell construction. 

On the other hand, Heiselberg et al. pointed out the importance of the inclusion 
into consideration of the finite size effects to realize the structured mixed phase [ 
They demonstrated that for rather large values of the surface tension the mixed phase 
constructed of rather small size droplets becomes energetically unfavorable compared 
to the two bulk-phase structure given by the Maxwell construction. Recently it was 
argued that the Gibbs condition of the equality of the electron chemical potentials is 
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not satisfied at the interface between tlie color-flavor-locked phase (no electrons) and the 
nucleon matter phase and it is described by the Maxwell construction [ |^ . In the above 
treatments the charged particle density profiles are uniform in each phase and furthermore 
corrections due to charge screening effects are disregarded. In spite of the question on 
the role of the screening has been raised long ago, up to now there was no consistent 
treatment of these inhomogeneity effects. Therefore, a further study of the screening and 
surface effects seems to be of prime importance. 



2. DENSITY FUNCTIONAL METHOD 

Consider a quark droplet (I) immersed in the (n, p, e) nuclear matter (II) within the 
Wigner-Seitz cell. The narrow boundary region between them is approximated by the 
sharp boundary dD. Then the total energy functional E is given by a functional of 
particle densities Pi[ H: 

E[p{r)] = I rfVeL+str(p[(r)) + / rfVe^J„+,,,(pf (r)) + / dhes + By, (1) 

where the first and second terms are the kinetic plus strong interaction contributions in 
both phases, the third term is the surface energy contribution. The last term is the usual 
Coulomb interaction energy, 

Ev = \j j dhd\' E (2) 

with Qi being the particle electric charge. Then the generating functional under the con- 
straint of the particle number conservation is introduced via the Legendre transformation, 

fi[p(r)] = E[p{v)] -T^f^lJ, d'rpi (r) - J] pf / ^Vpf (r). (3) 

i i 

Note that fi[p(r)] is the analog of the thermodynamic potential and the parameters p" are 
the chemical potentials. When we consider two conservation laws relevant in the mixed 
phase : baryon number and charge conservation, these quantities are well defined over 
the Wigner-Seitz cell, not restricted to each domain. Accordingly the baryon number 
and charge chemical potentials (pe and pg) , being linear combinations of p", become 
constants over the whole space, 

f^B = f^B = f^B, l^Q = /ig = /ig; (4) 

e.g. pq = Pg and Pb = fJ-n = 2pd + p^. 

Equations of motion are given by the variational principle: 

^-f^ = or pt = ^^|±^ - iQt/e)V-ir), a = {/, //}, (5) 

where V""(r) is the electric potential. 
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and satisfies the Poisson equation, 

^V' = AireY^Q^pl (7) 

i 

Note here that the electric potential can be shifted by an arbitrary constant, , and 
theory should be invariant under the redefinition of chemical potentials in a proper way; 
V — > V — /if — > fj,f + (Qf /e)V^°. One convenient choice is = that means that 
the electron density can be expressed as pe = (/ie — VY/ (3vr^). Another convenient choice 
is = —fie corresponding to pe = —V^/ (37r^), see [ |T|. 

If we expand the charge density in the r.h.s. of the Poisson equation with respect to 
around a reference value up to the second order, we arrive at the equation 

AV^" = i7iJ2 <5>i(^" = K") + ^IV" (8) 
i 

with V" = V"' — V^. The coefficient is the inverse of the Debye screening length, 

.dp\ 



(9) 

V"=VS' 



Electrons (quarks) can be treated as having the uniform charge density profile only if the 
droplet size and the inter-droplet distance are essentially smaller than the screening length 
^e,D {\,d)- Otherwise, as we shall demonstrate it below, the screening effects become 
essential and the charged density profiles are not uniform. Estimates show that the 
screening length for electrons Ae,D > 13 fm is longer than that for the quarks, Xq^D — 5 fm. 

Pressure in each phase is constant owing to the equation of motion (|^), and the pressure 
balance condition at the boundary dD reads 

pi = pii + p^^ (10) 

where Pg denotes the surface contribution. In the following we assume the spherical 
geometry and presence of the sharp boundary parameterized by the constant surface 
tension a. Then Eq. (|10|) is equivalent to the extremum condition, dVt/dR = 0, where R 
is the droplet radius. Eqs. (^), ( jlOD are manifestation of the Gibbs conditions. 

3. RESULTS AND DISCUSSION 

In Fig. 1 we show the effective energy of the Wigner-Seitz cell per droplet volume 
associated with the inhomogeneity of the electric field profile, 6oj = 3Q[p{r)]/ (AnR^), up 
to the second order in V"' as a function of the scaled droplet radius C, = R/\]j- The 
curves are presented for two values of the concentration of phase I, / = (R/Rw)^, Rw 
is the radius of the Wigner-Seitz cell. Solid curves are for / = 1/2 and dashed curves 
are for a tiny value / = 0.01. The "e.m." curves show the partial contribution of 
the (screened) electric field to the energy (oc /(VV")^c/^r) plus the surface energy. The 
difference between the non-labeled curves and the "e.m." curves shows the important role 
of the secondary (correlation) effects coming from the dependence of the bulk energies 
(ekin+str) on the inhomogeneous distribution of the electric field. The ratio of the screening 
lengths ao = X^/^d assumed to be 1. 
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We vary the value of the surface tension, /?i = 3(t/(A^/3o) with the energy scale parameter 
j3o, admitting its uncertainties. E.g., for typical values /ig — 170 MeV for the electron 
chemical potential, /i„ ~ 1020 MeV for the baryon (neutron) chemical potential, etc — 0.4 
for the strong coupling constant, and ~ 120 -t- 150 MeV for the strange quark mass, 
we estimate i3q ~ 1.6m^. Then, with the value a ~ 1.3m^ we obtain jSi ~ 0.7, whereas 
with 0" ~ 10 MeV/fm^ ~ 0.14m^ we would get Pi ~ 0.08. The Coulomb curves (labeled 
by "C") show the corresponding results obtained neglecting screening effects. We see 

1/3 

that the "C" curves demonstrate pronounced minima ai ^ = Pi ? while there is no 
minimum for Pi > Pic if one includes screening effects. According to our figure Pic — 0.6 
(would be Pic — 0.03 for "e.m." curves, when correlation energy effects are disregarded). 
Only for Pi < 0.01 the minima at the "C" curves differ not essentially from the minima 
at other curves. Hence our results clearly show that the screening effects are important 
for all realistic values of a and that the structured mixed phase is proved to be prohibited 
due to the kinetic instability of the droplets induced by the screening effects provided the 
surface tension is larger than a critical value. In absence of the mixed phase our charged 
distributions describe the boundary layer between two separated phases existing within 
the double-tangent (Maxwell) construction. Consideration of non-spherical droplets (rods 
and slabs) does not change our conclusions [ 0. 
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Figure 1. Contribution to the effective energy per droplet volume due to inhomogeneous 
charge distributions versus scaled droplet radius for the droplet concentration / = 1/2 
(solid curves) and a tiny concentration / = 0.01 (dashed curves), /Sq is the scale-parameter 
for the energy, and f3i is the parameter of surface tension (see the text). The "C" curves 
are calculated ignoring screening effects, the "e.m." curves include the electric field en- 
ergy (with screening effect included) plus surface energy, ignoring however the correlation 
effects. 



